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Abstract 

Measurements of Rb, the ratio of the bb cross-section to the qq cross-section in e + e~ 
collisions, are presented. The data were collected by the OPAL experiment at LEP 
at centre-of-mass energies between 182 GeV and 209 GeV. Lepton, lifetime and event- 
shape information is used to tag events containing b quarks with high efficiency. The 
data are compatible with the Standard Model expectation. The mean ratio of the eight 
measurements reported here to the Standard Model prediction is 1.055 ± 0.031 ± 0.037, 
where the first error is statistical and the second systematic. 
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1 Introduction 



The cross-section for bottom-quark pair production in e + e annihilation relative to the 
hadronic cross-section, 

_ <r(e + e~ — > 7/Z — > bb) 
b cr(e + e~ — > 7/Z — ► qq) ' 

is a sensitive probe of the Standard Model [1]. Measurements of R h have been made at 
the Z peak and at higher energies [2]. At the Z resonance, where fermion-pair production 
is dominated by Z decays, measurements of Rb provide a precise determination of the 
ratio of the Z — > bb partial width to the hadronic width 1 R® = T^/Phad- This quantity is 
of particular interest because of its unique sensitivity to electroweak radiative corrections; 
while sensitive to the top-quark mass, its dependence on other parameters, for example 
the Higgs boson mass and the strong coupling constant, is negligible [2]. Above the Z 
peak, the pure Z cross-section decreases and the contributions of photon exchange and 
7-Z interference become important. Possible new physics at a high energy scale might 
manifest itself as a deviation from the Standard Model prediction. 

In this paper, measurements of i?b at energies above the Z resonance are presented. 
The data were taken by the OPAL detector at the LEP e + e~ collider, at centre-of-mass 
energies, \/s, ranging from 182 GeV to 209 GeV, during the LEP2 programme. 

Above the Z peak a significant fraction of the observed fermion-pair events comes 
from radiative return to the Z through initial-state photon radiation. Only non-radiative 
events are considered here, according to the definition used by OPAL in the analysis of 
fermion-pair production at LEP2 [3]: 

• An effective centre-of-mass energy y/s' > 0.85^, where s' is defined as the square 
of the mass of the 7/Z propagator. 

1 R^ ) is the partial width ratio for Z decays and not the cross-section ratio measured in this paper. It 
is smaller than the cross-section ratio at the peak of the Z resonance by 0.0002. 
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• The predicted contribution of interference between initial- and final-state photon 
radiation is removed. 

At high energies, additional background sources also arise, mainly from W + W~ and 
ZZ decaying to four-fermion final states. 

Because of the limited statistics due to the low cross-section at LEP2, the double-tag 
technique as used at the Z resonance [4] is not optimal. Instead the measurement reported 
here relies on a single-tag method, and uses a sophisticated tagging algorithm based on 
lepton, lifetime and event-shape information to identify bb events. This algorithm is more 
efficient and has a higher purity than the one used in our previous measurement [5] , where 
i?b was measured up to 189 GeV. Hence the 183 GeV and 189 GeV results reported here 
supersede the previous values, while the measurements at higher energy are reported here 
by OPAL for the first time. 



2 The OPAL detector, data sample and simulation 

A detailed description of the OPAL detector can be found elsewhere [6]. For this analysis, 
the most relevant parts of the detector are the silicon micro-vertex detector, the track- 
ing chambers, the electromagnetic and hadronic calorimeters, and the muon chambers. 
The micro-vertex detector is essential for the reconstruction of secondary vertices. The 
central detector provides precise measurements of the momenta of charged particles from 
the curvature of their trajectories in a magnetic field of 0.435 T. In addition, it allows 
the identification of charged particles through a combination of the measurement of the 
specific energy loss dE/dx and the momentum. The electromagnetic calorimeter consists 
of 11 704 lead glass blocks, which completely cover the azimuthal 2 range up to polar an- 
gles of | cos#| < 0.98. Almost the entire detector is surrounded with four layers of muon 
chambers, after approximately one metre of iron from the magnet return yoke, which is in- 
strumented as a hadron calorimeter. Luminosity is determined using small-angle Bhabha 
scattering in the forward calorimeter [3]. 

2.1 Data sample 

The data used in this analysis were collected by the OPAL detector at LEP during 1997- 
2000, at centre-of-mass energies between 182 GeV and 209 GeV. They are grouped into 
eight energy points. The centre-of-mass energies and integrated luminosities of the data 
samples are summarized in Tabled 

2 The OPAL coordinate system is denned as a right-handed Cartesian coordinate system, with the x 
axis pointing in the plane of the LEP collider toward the centre of the ring, the z axis in the direction of 
the electron beam, and 9 and <fi defined as the usual spherical polar coordinates. 
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Nominal Energy 


Energy Range (GeV) 


(V5> (GeV) 


Integrated Luminosity (pb 1 ) 


183 


181.7- 184.0 


182.7 


56 


189 


188.5 - 189.1 


188.6 


185 


192 


191.4 - 191.7 


191.6 


29 


196 


195.4 - 195.8 


195.5 


77 


200 


199.3 - 199.9 


199.5 


78 


202 


201.5 - 202.0 


201.6 


36 


205 


204.5 - 205.5 


205.3 


74 


207 


205.5 - 208.6 


206.8 


137 



Table 1: The nominal energy, energy range, luminosity- weighted mean centre-of-mass 
energy and integrated luminosity of the data used in the analysis. 



2.2 Simulation 

Monte Carlo (MC) simulation is used to estimate the selection efficiency for bb events 
and to evaluate backgrounds. MC samples are generated at the nominal centre-of-mass 
energies given in Table 1. The KK2F program [7] is used to simulate e + e~ — > qq(7) events. 
Four-fermion background events are simulated with the grc4f generator [8] or with the 
KoralW [9] program with the grc4f matrix elements. The hadronization is performed by 
the Pythia 6.150 [10] string model, and heavy quark fragmentation is modelled according 
to the Peterson fragmentation scheme [11] with parameters tuned according to [12]. All 
the MC samples are passed through a detailed simulation of the OPAL detector [13]. 

The Standard Model cross-section predictions are calculated using Zfitter 6.36 [14]. 
Zfitter is used to calculate the Standard Model value of R c , the ratio of the cross-section 
for charm-quark pair production to the hadronic cross-section, at which the resulting 
is quoted, and to calculate i?b for comparison of the theoretical predictions with the 
measured values. Zfitter is also used for estimating the effects of interference between 
initial- and final-state photon radiation, as described in Section l3~4l 



3 Analysis procedure 

The measurement of i?b starts by selecting non-radiative hadronic events, applying fiducial 
cuts and rejecting four-fermion background events. The number of selected events, Ahad, 
is corrected for background and for the efficiency of the selection cuts to yield a corrected 
number of non-radiative hadronic events, A q q. The background includes the residual 
four-fermion background and feedthrough of events with lower effective centre-of-mass 
energy, \f^. The contribution from interference between initial- and final-state radiation 
is also removed. A b-tagging algorithm is applied to the selected events, and the number 
of tagged events, A tag , is similarly corrected for the same background contributions and 
interference to obtain A b ^, which still contains contributions from cc and light-quark-pair 
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events. is then calculated using: 

d N bb/ N <n ~ e c R c ~ (! ~~ ^c)e u ds n \ 

-Kb = , l-LJ 

— Cuds 

where eb, c ,uds are the efficiencies for non-radiative bb, cc and light-quark-pair events to 
pass the selection criteria and the b-tagging algorithm. 

In the following sections, we describe the selection of high quality non-radiative hadronic 
events, the rejection of four fermion background (Section 13. 1)1 . and the tagging algorithm 
used to select events originating from e + e~ — > bb ( Section 13. 2|) . This is followed by the 
description of the residual background estimation ( Section 13. 3J) . The correction for inter- 
ference is discussed in Section 13 A\ and finally the calculation of is presented (Section 
13.5)1 . The systematic uncertainties are described in Section 01 and the final results and 
conclusion are given in Section El 

3.1 Hadronic event selection 

Hadronic events, e + e~ — > j/Z — > qq, are selected based on the number of reconstructed 
charged particle tracks and the energy deposited in the calorimeters [15]. Further require- 
ments follow, demanding at least seven tracks and that the polar angle of the thrust axis, 
9 thr , satisfy | cos^^l < 0.9; the thrust axis is calculated from all tracks and clusters and 
corrected for double-counting as in [16]. The effective centre-of-mass energy, y/s', of the 
e + e~ collision is estimated as described in Ref. [3], and non-radiative hadronic events are 
selected if \J s' / s > 0.85. For the centre-of-mass energies analysed in this paper, the pro- 
duction cross-section of W + W~ events is comparable to that of non-radiative qq events. 
The production cross-section of ZZ events is about an order of magnitude smaller. Using 
the same techniques as in [17], events are rejected if they are identified as W + W~ events. 

The number of events at each energy selected by the above cuts (iVhad) is given in 
Table El Typical efficiencies, including the ss 90% efficiency of the cos# t hr cut, are 77% 
for bb and 80% for lighter (u, d, s and c) qq events. About 14% of the iV had sample are 
expected to be bb events. The selection of non-radiative events is less efficient for bb final 
states than for other flavours because of the generally larger missing energy due to the 
neutrinos in semileptonic b and c hadron decays. The remaining four-fermion background 
is estimated to be about 6% of the selected hadronic events. 



3.2 b-tagging 

In e + e~ — > 7/Z — > bb events, the quark and the anti-quark are typically boosted in 
opposite directions, and the subsequent hadronization is largely independent. In this 
analysis, each event is divided into two hemispheres defined by a plane that is orthogonal 
to the thrust direction. 

We use a hemisphere-based b-tagger designed for LEP2 Higgs searches [18], where the 
tagging of bb events is based on three nearly independent properties of the b-hadron and 
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its decay products: a lepton from a semileptonic b-hadron decay, the long lifetime, and 
kinematic differences between b-hadron decays and fragmentation in uu, dd, ss events, due 
to the hard fragmentation of the b-hadron and the high multiplicity of its decay products. 

For the lepton tag, semileptonic b-hadron decays are identified from their resultant 
electron or muon and the lepton momentum is used as the b-tag variable. Electrons are 
identified using an artificial neural network (ANN) [4] while muons are identified using 
information from the muon chambers in association with the tracking chambers, as in [19]. 
Electrons from photon conversions are rejected using an additional ANN described in [20]. 

For the lifetime tag, secondary vertices are reconstructed. Variables based on the 
vertex significance (the distance between the reconstructed primary and secondary ver- 
tices divided by the corresponding uncertainty) and the impact parameters of the tracks 
associated with the vertex are combined using an ANN scheme as in [18]. 

For the kinematic tag, three variables are combined using another ANN: the number of 
detected particles in the central part of the hemisphere, the angle between the hemisphere 
axis and its boosted sphericity axis, and the C parameter [21], the latter two being 
calculated in the rest frame of all the particles associated with the hemisphere. 

In each hemisphere, the output of the lepton tag, the lifetime ANN and the kine- 
matic ANN are combined with an unbinned likelihood calculation [22] , and the likelihoods 
£hemii5 ^ncmii' ^hemii' ^hemi2 5 ^Lmi25 ^hemi2 are obtained. The likelihood C\ erail is the prob- 
ability for the first hemisphere to contain a b quark. C\ eral2 ls the same probability for 
the second hemisphere, and £himii(2) are ^ ne corresponding probabilities for c-quark (light 
quarks). The two hemispheres' outputs are then combined into a single event b-tagging 
likelihood variable £ e vent, using 

r , rb rb 
r> ' b ■ < ~hemil' ! -'hcmi2 



event r , rb rb _i_ „ re rc < r , ruds ruds ' 

' b / -hemil A -'heini2 ' ' c *-hemil' i -"hemi2 ' ' uds *-hemil" i -hemi2 

The normalization parameters are set to rb = 0.165, r c = 0.253 and r u d s = 0.582, which 
are the Zfitter theoretical predictions for R h , R c and R u d s at about 190 GeV. The 
same values are used for all data sets. The number of events, iV ta g, satisfying a b-tagging 
cut of £ e vent > 0.3 is determined; this cut value minimizes the total uncertainty in the 
measurement of Rb- Typical efficiencies for the b-tagger in the hadronic event sample are 
65%, 6.3% and 1.5% for bb, cc and light-quark pair events respectively; for the residual 
four-fermion events the tagging efficiency is about 5%. On average, 77% of the tagged 
sample is expected to be bb events. The distribution of £ ev ent for the 182-209 GeV data is 
shown in Figure ^ together with the expectation from the Monte Carlo. Good agreement 
between the data and the Monte Carlo is observed. 



3.3 Background estimation 

The residual four-fermion background is estimated from Monte Carlo simulation and 
subtracted from the number of selected hadronic (iV had ) and b-tagged (-/V tag ) events. It 
is 6.5% for the hadronic event sample and 2.5% for the b-tagged sample. The dominant 
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Figure 1: £ CV cnt distributions for selected non-radiative hadronic events. The points rep- 
resent the data (182-209 GeV) and the histogram represents the Monte Carlo prediction, 
normalized according to the integrated luminosity of the data, with the open area rep- 
resenting the expected bb content and the hatched areas representing the contributions 
from cc light-quark pair and four-fermion events. The arrow indicates the cut used. 
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contribution to the four-fermion background originates from W + W and ZZ final states, 
with only 5% coming from non-W + W~ and non-ZZ events. 

Backgrounds from two-photon collisions and e + e~ — > t + t~ events were found to be 
negligible. The selected event sample has a small contamination from radiative hadronic 
events where \J s' /s was overestimated. This contamination is estimated using Monte 
Carlo simulation and subtracted from A^ a d and iVt ag when calculating iVqq and iV bb -. 
Since both the b-tagged and the hadronic samples include such contamination, the effect 
on the measured value of i?b is small. About 3% of the selected hadronic events and 4% 
of the b-tagged events have a true value of \/lP smaller than 0.85^. 

The contributions from non-radiative cc and light-quark-pair events that are b-tagged 
are taken into account by the two negative terms in the numerator of equation 1. The 
secondary production of heavy quarks, commonly referred to as gluon splitting, can also 
lead to primary light-quark-pair events being selected by the b-tagging algorithm. This 
contribution is fully accounted for by the efficiencies in equation 1 if the rates of this 
production in the Monte Carlo are correct. These rates were measured at LEP1 [23-25]. 
We extrapolate the measured values to the higher energies using theoretical predictions 
given in [10,26,27]. We obtain production rates per hadronic event, g C c = 0.064±0.010 and 
Sbb = 0.0065±0.0015 at y/s = 200 GeV where the uncertainties include the uncertainty on 
the measurements and the variation of the theoretical models used for the extrapolation 
to higher energy. The default Monte Carlo values are g C c = 0.032 and g bb = 0.0053. We 
correct the value of g cC in the Monte Carlo according to the extrapolated value. Since 
the Monte Carlo value of g bb agrees with the data we do not correct the Monte Carlo for 
g bb but take the difference between the Monte Carlo value and the maximal extrapolated 
rate (0.0053 — 0.0080) as the corresponding systematic uncertainty. 

3.4 Subtraction of interference between initial- and final-state 
photon radiation 

The interference cross-section is defined as the difference between the cross-section in- 
cluding interference between initial- and final-state photon radiation and the cross-section 
excluding interference. It is calculated using Zfitter with the Standard Model value of 
the forward backward asymmetry, for e + e~ — > bb, e + e _ — > cc and e + e~ — > uu, dd, ss, in 
bins of cos 9 where 9 is the angle between the fermion direction and the direction of the e~ 
beam. The interference cross-section can be positive or negative, depending on the value 
of 9. The thrust axis gives a good estimate of the fermion direction, so we assume that 
the interference correction estimated in bins of cos# can be applied in bins of cos#thr- The 
fractions of b, c and uds events in the hadronic and b-tagged samples are estimated from 
Monte Carlo samples generated excluding interference. The interference cross-section for 
each flavour is then weighted according to the estimated flavour composition, and the 
resulting numbers of events are subtracted from the number of b-tagged events and from 
the number of qq events to obtain iV bb and N qq , which appear in equation 1. These cor- 
rections are typically 0.3% for the hadronic sample and —0.5% for the b-tagged sample. 
The uncertainty on the above procedure is discussed in Section 14. 1.31 and given in Table 01 
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Data set 


iVhad 


iV tag 




^bb 


£b 


e c 


^uds 


i? c SM 


i?b ± stat 


183 


1172 


185 


1199.8 


176.1 


0.60 


0.062 


0.012 


0.2532 


0.207 ±0.018 


189 


3496 


447 


3561.1 


422.7 


0.58 


0.059 


0.012 


0.2551 


0.165 ±0.010 


192 


573 


80 


584.7 


76.4 


0.58 


0.064 


0.014 


0.2560 


0.174 ±0.025 


196 


1343 


182 


1370.6 


172.7 


0.56 


0.060 


0.014 


0.2570 


0.181 ±0.017 


200 


1312 


166 


1333.9 


157.3 


0.58 


0.061 


0.014 


0.2581 


0.164 ±0.016 


202 


591 


72 


604.6 


68.3 


0.56 


0.062 


0.014 


0.2585 


0.154 ±0.024 


205 


1124 


136 


1140.6 


128.0 


0.55 


0.061 


0.014 


0.2592 


0.158 ±0.018 


207 


2046 


251 


2078.5 


237.5 


0.54 


0.057 


0.014 


0.2597 


0.169 ±0.014 



Table 2: Number of selected (iV had , iV tag ) and corrected (iV q q, N h ^) events; average (over 
all Icos^thrl < 0.9 bins) selection (including b-tagging) efficiencies; the Standard Model 
prediction for R c and the resulting value of (with statistical uncertainties only) for the 
eight energies. Note that R^ is calculated from the values obtained in bins of cos 9 t hr and 
is not a direct calculation using the numbers in the table. 



3.5 Calculation of 

To avoid dependence on the Monte Carlo angular distribution, the analysis is performed 
in bins of Icos^thrl- In each |cos#thr| bin, the efficiencies and corrected numbers of qq 
and bb events are determined. The efficiencies, eb, c ,uds, are defined as the number of MC 
events of the specific flavour b-tagged in the | cos (9thr| bin, divided by the number of events 
of that specific flavour that were generated in that bin. 

_Rb is then calculated in each bin in the range |cos# t hr| < 0.9. Figure El shows the 
value of i?b versus | cos^thrl for the 189 GeV and 207 GeV samples. The final value of R^ 
is determined by taking the average over all bins. The total numbers of selected hadronic 
and tagged events, the values of the variables in equation 1, as well as the resulting R^ 
value with its statistical uncertainty, are summarized in Table El 

The resulting i?b values are quoted for the case where R c is set to its Standard Model 
value, R^ M , given in Table El The dependence of the result on the assumed value of R c 
can be parametrized as: 

AR h = b (R c - R™) . 

The parameter b, derived from equation 1, is determined separately for each centre-of-mass 
energy and is given in Table 4. 



4 Systematic uncertainties 

The systematic uncertainties considered for these measurements are described in the fol- 
lowing paragraphs. A breakdown of the systematic uncertainty at 207 GeV is given in 
Table El and the total systematic error is given for all energies in Table 4. 
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Figure 2: Measured values of i?b for: a) the 189 GeV sample; b) the 207 GeV sample 
with statistical uncertainties only. The dashed line is the average value and the solid lines 
show the the band corresponding to the one standard deviation uncertainty. 
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Uncertainty source 


AR h /R h (%) 


b fragmentation 


0.6 


b decay multiplicity 


0.6 


b hadron composition 


0.3 


b lifetime 


0.1 


c fragmentation 


0.8 


c decay multiplicity 


0.1 


c hadron composition 


0.3 


c lifetime 


0.1 


semileptonic branching ratio 


0.1 


semileptonic decay model 


0.1 


Total b,c physics modelling 


1.3 


K°, A rate 


0.2 


Interference 


0.3 


Four-fermion background 


0.1 


gbb 


0.1 


gcc 


0.1 


Total physics modelling 


1.3 


Track reconstruction 


3.2 


Lepton ID 


0.5 


Non-radiative event selection 


0.1 


Monte Carlo statistics 


0.6 


Total detector effects 


3.3 


Total systematic uncertainty 


3.5 



Table 3: Systematic uncertainty breakdown at y/s = 207 GeV. The uncertainties at other 
centre-of-mass energies are similar. 

4.1 Physics modelling 

4.1.1 Bottom and charm physics modelling 

Uncertainties in bottom and charm fragmentation and decay properties are discussed 
below. The variation of parameters is realized by reweighting Monte Carlo events to the 
modified distribution. 

• b fragmentation: Although the mean scaled energy 2{E\>)/y/s of weakly decaying 
b hadrons is expected to change from LEP1 to LEP2 energies, the free parameter 
e° P of the Peterson fragmentation function [11] used in the default MC for heavy 
flavour events is assumed not to vary with energy. The value in bb events is varied 
in the range 0.0030 < ep < 0.0048, which corresponds to a variation of the mean 
scaled energy 2(E\>)/y/s of weakly decaying b hadrons in Z decays in the range of 
2(Eb)/y/s = 0.702 ± 0.008 [12]. In addition, the heavy-quark fragmentation model 
is changed to that suggested by Collins and Spiller [28] and to that of Kartvelishvili 
[29]. The largest difference in all the above tests is from the Peterson parameter 
variation and is taken as the systematic uncertainty. 
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• b decay multiplicity: The multiplicity of the charged decay products of hadrons 
containing a b quark is varied in the Monte Carlo simulation according to [12]. 

• b hadron composition: The tagging efficiency differs for the various b hadron 
species. The fractions of b hadrons and their uncertainties have been taken from [30]. 
The fractions /(B°) and /(bbaryon) are varied independently within ranges of the 
experimental measurements, and their variation is compensated by the fraction of 
B° and B^ mesons. 

• b hadron lifetimes: The lifetimes of the different b hadrons are varied in the 
Monte Carlo by their uncertainties according to [30]. 

• c fragmentation: Simulated cc events have their Peterson fragmentation parame- 
ter e c p varied in the range of 0.022 < e c P < 0.039. This corresponds to a variation of 
the mean scaled energy 2(E c )/y/s of weakly decaying c hadrons in Z decays in the 
range of 2{E c )/y/s = 0.484 ± 0.010 [12]. 

• c decay multiplicity: The average charged particle multiplicities of D + , D° and 
D+ decays are varied in the Monte Carlo within the ranges of the experimental 
measurements [31]. 

• c hadron composition: The D° fraction is written as /(D°) = 1 — /(D + ) — /(D+) — 
/(cbaryon)- The last three parameters are varied independently by their uncertainties 
according to [32] to evaluate the uncertainty on the charm efficiency. 

• c hadron lifetimes: Charmed hadron lifetimes are varied within their experimental 
uncertainties according to [30]. 

• semileptonic decay modelling and branching ratios: The semileptonic branch- 
ing ratios of bottom- and charm-hadrons are varied within their experimental un- 
certainties [12]. The uncertainty due to the semileptonic decay model used in the 
simulation is estimated as in [24]. 

4.1.2 Inclusive K° and A production rate 

The total production rates of K° and A in the Monte Carlo are varied by 4% and 7% 
respectively according to [4]. This variation contributes an uncertainty of 0.18%. 

4.1.3 Initial- final state interference 

To asses the uncertainty on the correction for interference, we used the results of a study 
made in [3] by replacing the Zfitter predictions with those of the KK2f program [7]. 
We assign the difference between the two predictions as a systematic uncertainty. 

4.1.4 Four-fermion background 

The uncertainties on the measured W + W~ and ZZ cross-sections [2] and the uncertainty 
on the luminosity are found to have a negligible systematic effect on R h . The background 
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from W- and Z-pairs has the highest probability to be accepted in the tagged sample when 
one or both bosons decay into a final state containing a charm- or bottom-quark. The 
systematic uncertainty on the W- and Z-pair tagging efficiency is estimated by varying 
the tagging efficiency by 5%. This variation introduces an uncertainty of 0.15%. The 
effect of the detector resolution is also taken into account, as described in Section ^. 2. 11 

4.1.5 Secondary heavy-quark production 

As discussed in Section [%!B1 the difference between the Monte Carlo rate for secondary 
bb production and the highest expected value is taken as an uncertainty. Typically, this 
process contributes a relative uncertainty of about 0.1%. The secondary cc production 
rate is taken from [23,24] and weighted according to [26] as described in Section l3~3l The 
uncertainty on this estimation contributes a relative uncertainty of 0.13%. 

4.1.6 Monte Carlo bb production 

The transition iVhad - ► N qq has a small dependence on the Monte Carlo value of R\>. The 
default Monte Carlo value of i?b is varied within the measured uncertainty (Table 2). The 
effect on R h is negligible. 

4.2 Detector effects 
4.2.1 Track reconstruction 

The detector performance in the Monte Carlo is varied as in [4, 18]. The effect of the 
detector resolution on the track parameters is estimated by degrading the resolution of 
all tracks in the Monte Carlo simulation. This is done by applying a single multiplicative 
scale factor to the difference between the reconstructed and true track parameters. A 5% 
variation is applied independently to the rip and rz track parameters. These variations 
together contribute an uncertainty of 1.8%. In addition, the efficiency for assigning mea- 
surement points in the silicon microvertex detector to the tracks is varied by 1% in r<p 
and 3% in rz. Each of these variations contributes an uncertainty of 1.3%. The track 
reconstruction efficiency is varied by 1%. This is done by randomly discarding 1% of the 
tracks, and contributes an uncertainty of 1.8%. The alignment of the silicon microvertex 
detector is changed in the Monte Carlo as in [4]. This variation gives an uncertainty of 
0.4%. The systematic uncertainties resulting from the individual variations are summed 
in quadrature. 

Since the largest part of this uncertainty is due to the effect on e u d s , we verify the 
validity of the 5% variation in the rip track parameters by using events that are rejected 
by the b-tagger. We look at the 182-209 GeV sample as well as at data collected in the 
same years at the Z peak for detector calibration. Events with a low value of £hemi in 
one hemisphere are selected, giving a sample dominated by light-quark pairs. The value 
of the vertex significance in the opposite hemisphere is then examined. Figure El shows 
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the fraction of events with vertex significance less than a negative cut value in the data 
and in the Monte Carlo with 5% variation in the r<p track parameters, normalised to the 
fraction in Monte Carlo. This study suggests that the difference between the data and 
the Monte Carlo is smaller than the variation on the Monte Carlo and we conclude that 
the 5% variation is adequate. 
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Figure 3: Ratios of fractions of hemispheres with vertex significance smaller than the cut 
value on the abscissa. The data (open squares) and the Monte Carlo with 5% deterioration 
of the tracking parameters resolution in rep (solid triangles) are normalised to the default 
Monte Carlo. The uncertainties are a combination of the data and Monte Carlo statistical 
uncertainties. Hemispheres are selected if the opposite hemisphere in the event has a low 
value of £hemi, giving a sample 75% pure in uds events. The 182-209 GeV sample is 
shown in a) and the 91 GeV sample in b). 
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4.2.2 Lepton ID 



The electron identification efficiency is varied by 15% and the electron fake rate by 30%. 
These values are 50% larger than those used in [5] following a later comparison of OPAL's 
data and Monte Carlo. The muon identification efficiency is varied by 2% and the fake 
rate by 14%. This source contributes typically an uncertainty of 0.45% on i?b which is 
predominantly from the variation on the electron identification efficiency. 

4.2.3 Non-radiative event selection 

We repeated the analysis using a different method for estimating s', described in [33]. 
The resulting i?b was found to be different by less than 0.1%. This comparison is almost 
free of statistical uncertainties; thus we take this difference as an uncertainty. 

4.2.4 Monte Carlo statistics 

The finite size of the Monte Carlo samples contributes between 0.6% and 0.8% uncertainty 
to the systematic error. 

4.3 Cross-checks 

As a cross-check, the analysis is repeated on data collected in the same years at the Z peak 
for detector calibration. A value of R\,(y/s = mz) = 0.210 ± 0.001(stat.) ± 0.006(syst.) is 
obtained. The systematic uncertainty is determined as for the high-energy data samples. 
This result is consistent with the LEP1 average value of = 0.21643 ± 0.00073 [2]. 

Another cross-check is made by comparing the likelihood of the hemisphere to contain 
a b-quark in Monte Carlo and data. This comparison is done for both the 182-209 GeV 
data and the calibration data at the Z peak. We select hemispheres where the b-tagger 
output in one hemisphere indicates a high probability for the event to be a bb event, and 
look at the opposite hemisphere. The samples are 60% and 88% pure in bb events for the 
high energy and the Z peak data respectively. We also look at hemispheres in events that 
are rejected by the b-tagger in the opposite hemisphere, giving samples with only 6% bb 
events. The results of this cross-check for the 182-209 GeV data and for the calibration 
data at the Z peak are shown in Figure EJ The small excess of data with respect to the 
Monte Carlo for £homi > 0.85 in the anti-b-tagged plot of the 182-209 GeV sample (Figure 
IU (b)) is covered by the detector resolution uncertainty. It is consistent with the result 
of the study shown in Figure El (a). The difference between the Monte Carlo and the 
data in the b-tagged 91 GeV plot (Figure |H (c)) is also covered by the detector resolution 
uncertainty. 
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Figure 4: b-tagger output from one hemisphere for events selected by the opposite hemi- 
sphere. The top plots are for all data events in the range 182-209 GeV, and the bottom 
two are for the Z sample, a) and c) are for b-tagged opposite hemispheres and b) and d) 
are for anti-b-tagged opposite hemispheres. 



5 Results and conclusion 



The e + e~ — > bb production rate per non-radiative e + e~ — > qq event has been measured 
using data collected by the OPAL detector at LEP at centre-of-mass energies between 182 
GeV and 209 GeV. The results are summarized in Table 4, together with earlier OPAL 
results from 130-172 GeV [5]. The measurements are shown in Figure where they 
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are compared with the predictions of the Standard Model, calculated using Zfitter. 
Good agreement between the Standard Model and the measurements is observed. A 
comparison of the eight measurements reported here to the Standard Model prediction 
gives: x 2 — 5-0, or 76% probability of obtaining a larger difference than observed. Since 
the x 2 test only uses the absolute difference between the prediction and measurement, we 
also made another test assuming that the ratio of the measurement to the prediction is 
constant. We perform a x 2 fit f° r t ne ratio using only the statistical uncertainties and 
add the systematic uncertainties, which are assumed to be fully correlated, to the result. 
We obtain a ratio of 1.055 ± 0.031 (stat.) ± 0.037 (sys.). Both tests suggest that these 
measurements are consistent with the Standard Model. 

The measurement technique presented here has a larger selection efficiency and purity 
and thus yields a smaller statistical uncertainty and a smaller dependence on R c than the 
one used previously. The 183 GeV and 189 GeV results are in agreement with [5] and 
supersede the previous measurements. 
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Figure 5: The measured R h values (points with error bars) and the Standard Model 
prediction (solid line). 
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(GeV) 


R h ± stat ± sys 


h — Ail * 

u {Rc-R™) 


133.3 [5] 


0.190 ±0.023 ±0.007 


-0.12 


161.3 [5] 


0.195 ±0.035 ±0.007 


-0.11 


172.1 [5] 


0.091 ±0.034 ±0.005 


-0.11 


182.7 


0.207 ±0.018 ±0.007 


-0.085 


188.6 


0.165 ± 0.010 ± 0.006 


-0.082 


191.6 


0.174 ± 0.025 ± 0.006 


-0.088 


195.5 


0.181 ±0.017 ±0.006 


-0.092 


199.5 


0.164 ± 0.016 ± 0.006 


-0.083 


201.6 


0.154 ± 0.024 ± 0.005 


-0.088 


205.3 


0.158 ± 0.018 ± 0.006 


-0.088 


206.8 


0.169 ± 0.014 ± 0.006 


-0.082 



Table 4: i?b values with statistical (first) and systematic uncertainties. The last column 
gives the dependence of the i?b value on R c . No variation of R c is included in the sys- 
tematic uncertainty in the second column. The first three results are from the previous 
OPAL measurement [5]. 
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